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ABSTRACT 

We discuss near-IR images of six passive galaxies (SSFR< IIP 2 Gyr -1 ) at redshift 1.3 < 2 < 

2.4 with stellar mass M ~ 10 11 M 0 , selected from the Great Observatories Origins Deep Survey 
(GOODS), obtained with WFC3/IR and the Hubble Space Telescope ( HST ). These WFC3 images 
provide the deepest and highest angular resolution view of the optical rest-frame morphology of such 
systems to date. We find that the light profile of these, galaxies is generally regular and well described 
by a Sersic model with index typical of today’s spheroids. We confirm the existence of compact 
and massive early-type galaxies at 2 ~ 2: four out of six galaxies have r e ~ 1 kpc or less. The 
WFC3 images achieve limiting surface brightness //. ~26.5 mag arcseP 2 in the F160W bandpass; 
yet there is no evidence of a faint halo in the five compact galaxies of our sample, nor is a halo 
observed in their stacked image. We also find very weak “morphological k-correction” in the galaxies 
between the rest-frame UV (from the ACS 2 band), and the rest-frame optical (WFC3 H band): the 
visual classification, Sersic indices and physical sizes of these galaxies are independent or only mildly 
dependent on the wavelength, within the errors. 

Subject headings: cosmology: observations — galaxies: fundamental parameters — galaxies: evolution 


1. INTRODUCTION 

Most of the stellar mass observed in today’s early-type 
galaxies has assembled at high redshifts (z > 2, see Ren- 
zini 2006). Thus, the evolution of the mass functions 
and morphology of these galaxies is expected to provide 
information on the mechanisms that have assembled and 
shaped them in their current state, e.g. merging and 
accretion. 

Recent works agree that the most massive early-type 
galaxies (with M/M© > 10 11 ) were already in place 
at z~ 1, and characterized by a mass-dependent evo- 
lution of the luminosity function, whereas galaxies with 
lower mass evolve faster (Thomas et al. 2005; Frances- 
chini et al. 2006, Bundy et al. 2006, Scarlata et al. 2007). 

Galaxies at 1.5 < 2 < 3 with stellar mass and 
SEDs similar to those of local early types have also 
been identified and studied thanks to a number of ded- 
icated color-selection criteria (e.g. Franx et al. 2003; 
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Daddi et al. 2004; Kong et al. 2006; Ubert et al. 2009; 
McCracken et al. 2009), as well as spectroscopic sur- 
veys (e.g. Cimatti et al. 2008; Kriek et al. 2008). 
These galaxies have photometrically-measured stellar 
mass m ~ 10 11 M© and specific star formation rate 
SSFR< 10~ 2 Gyr -1 , nearly one order of magnitude lower 
than that of the Milky Way. There is broad agree- 
ment that the number density of these high-redshift 
“elliptical galaxies” rapidly increases from z~ 2 to 
z~ 1 (Fontana et al. 2006; Arnouts et al. 2007; II- 
bert et al. 2009). The mechanisms through which they 
build up their mass, however, are not yet known. Espe- 
cially intriguing is the recent discovery that at least some 
massive elliptical galaxies at 2 > 1.5 are, on average, 
~ 3x smaller, and thus « 30x denser, than their local 
counterparts (Daddi et al. 2005; Trujillo et al. 2006; Tru- 
jillo et al. 2007; Zirm et al. 2007; Longhetti et al. 2007; 
Cimatti et al. 2008, Van Dokkum et. al. 2008, 
Buitrago et al. 2008), a fact lacks an explanation 
in terms of evolutionary mechanisms (e.g. Hop- 
kins et al. 2009a). Recent works (Valentinuzzi et al. 2009; 
Trujillo et al. 2009), however, have identified similarly 
dense massive galaxies in the local universe, although 
their spatial density seems smaller than the 2 ~ 2 galax- 
ies. 

Some have suggested that current observations at high 
redshift could have missed low surface-brightness halos 
surrounding these galaxies, due to the relatively limited 
sensitivity of current imaging data (Hopkins et al. 2009b; 
Mancini et al. 2009). If true, this could help explain 
the apparently dramatic morphological transformation 
of passive galaxies from 2 ~ 2 to the present. Another 
possibility is that current surveys of the local universe 
(e.g. the SDSS, Stoughton et al. 2002) could also have 
missed a relatively large fraction of very compact, early- 
type galaxies due to the limited ground-based angular 
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resolution (1 arcsec corresponds to a physical scale of 
~ 1 kpc at 2 = 0.05). 

A crucial step is to obtain an unbiased and sta- 
tistically accurate measures of the size distribution of 
early-type galaxies at both high and low redshift (Cas- 
sata et al. in prep.). So far, while small high-redshift 
samples have been imaged with NICMOS at rest-frame 
optical wavelengths (Trujillo et al. 2007; Zirm et al. 2007; 
Longhetti et al. 2007; Van Dokkum et. al. 2008), most 
of the galaxies only have rest-frame UV images from 
ACS and the dependence of the morphology of these 
galaxies on the wavelength has not been characterized 
(Cimatti et al. 2008; Daddi et al. 2005). 

In this letter we use the unprecedented sensitivity and 
angular resolution of WFC3/IR images recently acquired 
in the HUDF to study the optical rest-frame morphology, 
as well as the spectral energy distributions, of a small but 
well defined sample of low SSFR galaxies at z~2. The 
WFC3 images, in the Y, J and H bands (F105W, F125W 
and F160W, respectively) have PSF with FWHM=0.18” 
at 16000A) and reach 1 — a surface brightness fluctua- 
tions of 27.2, 26.6 and 26.3 AB/arcsec 2 , respectively, the 
sharpest and deepest ever obtained to date. 

2. DATA ANALYSIS 

The imaging data comprise the first epoch of 
ultra-deep near-IR WFC3/IR observations acquired by 
Bouwens et al. (2009, program ID=G011563, P.I. G. 
Illingworth) in the HUDF (Beckwith et al. 2006), a 
« 3.3 x 3.3 arcmin 2 field fully embedded in the GOODS- 
South field (Giavalisco et al. 2004) and imaged with the 
same photometric system, i.e. the ACS BViz. This first 
epoch of data covers an area roughly equal to the foot- 
print. of the WFC3/IR camera (2.1 arcmin 2 ) and consists 
of images in the F105 (Y), F125W (J) and F160W (H) 
filters. We have carried out our independent reduction 
of the raw data and after rejecting images affected by 
persistence in the J band, our final stacks reach 1 — a 
surface brightness fluctuations of 27.2, 26.6 and 26.3 in 
the three bands, respectively. The WFC3 images have 
been drizzled to a 0.06” pixel scale. 

We selected a sample of high redshift (z> 1.3) pas- 
sive (SSFR< l()~' 2 Gyr~ 1 ) and massive (M> 10 10 M/M S ) 
galaxies from the GOODS multi-band GUTFIT photo- 
metric catalog by Grogin et al. (in preparation). This 
catalog contains matched photometry for all the photo- 
metric bands used for the GOODS multi-facility program 
from the U to the 8 pm IRAC band (the MIPS 24 /<m 
photometry is not included). All our sample galaxies, 
however, have no detection at 24/xm down to a 1 — <r 
limit of 5//Jy. Fontana et al. (2009) have shown that it is 
possible to select a sample of passive galaxies at high-z 
combining MIPS 24pm to the optical-NIR spectral en- 
ergy distribution. We used both the BzK color-selection 
criteria for passive galaxies by Daddi et al. (2004), as well 
as multi-band SED fitting to identify the sample. While 
the BzK technique provides a well characterized crite- 
rion to select galaxies with low specific star-formation 
rate (dubbed “pBzK” . or passive BzK), the method nat- 
urally has finite completeness and contamination, since 
photometric errors and the intrinsic dispersion of galax- 
ies' broad-band UV/ Optical SED scatter galaxies out 
of the selection window and interlopers into it. Thus, 
given the overall limited size of the sample, we have aug- 


TABLE 1 

The sample of passive galaxies 


ID 

Z W 

log(M/M 0 ) 

log(SSFR) 

[Gyr- 1 ] 

Sersic 

H-band 

r e 

[kpc] 

19389 

1.307p 

10.41 

-5.55 

2.99 ^ > 

1.02 W 

22704 

1.384s 

10.70 

-5.55 

2.65 

0.50 

23495 

2.349p 

11.14 

-3.39 

(3) 

< 0.38 

23555 

1.921s 

10.82 

-2.98 

1.97 

0.44 

24279 

1.980s 

10.63 

-3.39 

1.63 

0.37 

2-4626 

1.317s 

11.10 

-2.11 

7.42 

3.69 


Note. — (1): s and p indicate respectively spectroscopic and 
photometric redshifts; (2) the morphology of this object is repro- 
duced by adding to the Sersic model a central psf component; (3) 
the morphology of this object is reproduced by a single central psf 
component, so the size reported here is just an upper limit. 

mented it with galaxies selected for being massive and 
passive according to the above criteria based on fitting 
their observed broad-band GUTFIT SED to stellar pop- 
ulation synthesis models. We used the Chariot & Bruzual 
(2007) models with Salpeter IMF and lower and upper 
mass limits 0.1 and 100 M 3 , respectively, together with 
the Calzetti et al. (2000) obscuration law to account for 
the possible (modest) presence of dust and the Madau 
(1995) cosmic opacity. The results of the SED fitting are 
reported in Table 1, and the SEDs are shown in Fig. 1. 

The final sample, summarized in Table 1, consists of six 
galaxies, four of which have spectroscopic redshifts, i.e. 
22704 and 23555 from Cimatti et al. (2008), 24279 from 
Daddi et al. (2005) and 24626 from Vanzella et al. (2008), 
and two for which we derived accurate (Az/(1 + z) ~ 
0.05) photometric redshifts from the GUTFIT multi- 
band photometry. Four out of six galaxies fulfill the 
pBzK criterion. 

We have used GALFIT (Peng et al. 2002) to fit the 
light profile of the galaxies in the z, Y, J and H bands to 
the Sersic model 


7(r) = I e exp 





(1) 


where 7(r) is the surface brightness measured at dis- 
tance r, I e is the surface brightness at the effective ra- 
dius r e and b„ is a parameter related to the Sersic in- 
dex n. For n=l and n=4 the Sersic profile reduces re- 
spectively to an exponential or deVaucouleurs profile. 
Bulge dominated objects typically have high n values 
(e.g. n > 2) and disk dominated objects having n around 
unity. Cassata et al. (2005), Ravindranath et al. (2006) 
and Cimatti et al. (2008) showed that this method gives 
basically unbiased estimates of the Sersic index and ef- 
fective radius. 

The star used by GALFIT to convolve the model has 
been obtained averaging six bright unsaturated stars in 
the UDF/WFC3 field. GALFIT has been run consider- 
ing the sky as a free parameter and using a fitting region 
around each galaxy of 6 x 6 arcsec 2 . However, we have 
run tests using a different region size and fixing the sky 
to different values, and the Sersic indices and sizes so ob- 
tained were in agreement with the fiducial values within 
20 %. 


3. THE MORPHOLOGICAL PROPERTIES OF GALAXIES 
ATZ-2. 

Figure 1 shows for each galaxy the z— and 77— band 
images, the H-band residuals after subtracting the best 
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Fig. 1. — For each galaxy in the sample we show the postage stamp in z-band (1st column), in if-band (2nd column) and the residual 
after subtracting the best fit model in FT— band (3rd column). Each postage stamp is 6x6 arcsec 2 . In the 4th column, we show the galaxy 
profile in ii-band, overplotted to the profile of the PSF (dotted line) and a de Vaucouleurs profile with r e = 3 kpc (dashed line), together 
with the residual between galaxy and model profiles (bottom). Finally, in the 5th column we report the SED with the photometric best-fit 
model. 


GALFIT model and the model and galaxy profiles, to- 
gether with the SEDs. Recall that at z ~ 2, the z band 
probes the rest-frame UV and the H band rest-frame op- 
tical. 

Even on visual inspection, all the galaxies except 24646 
appear to be very compact. In Fig. 1 the Lf-band im- 
ages seem to be larger than the z-band ones, but this 
is just the result of the broader PSF (FWHM=0.16” in 
the H band vs. FWHM=0.12” in z). The light pro- 
files in both passbands are actually rather similar. At 


all wavelengths, all galaxies appear relatively symmetric, 
in close resemblance to present-day spheroidal galaxies. 
The only exception is galaxy 24626, the most extended 
of our sample, whose isophotes deviate from an ellipti- 
cal one. A low surface brightness shell-like structure is 
visible in the z-band image, on the lower-left side, while 
the position angle is different for different isophote lev- 
els. These irregularities can be the result of ongoing or 
recent interactions. 

The GALFIT Sersic indices and sizes from the H-band 
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images are listed in Table 1; they quantitatively validate 
the conclusions derived from the visual inspection. Of 
the six galaxies, four have Sersic index >2 (19389, 22704, 
23555 and 24626). In the case of object 19389 we had to 
add a central unresolved (PSF) component to the Sersic 
model to obtain a reasonable fit. This central stellar ob- 
ject contains less than 10% of the total light in the object, 
and may be indicative of an AGN. Galaxy 23495 also has 
a peculiar morphology: based on the FWHM, the object 
is barely resolved in all bands, Fig. 1 shows that the 
inner part of the galaxy’s profile (r< 1 arcsec) follows 
quite well the PSF profile, while the GALFIT best fit is 
that of an unresolved (PSF) light profile, again suggest- 
ing the presence of an AGN. Object 23495 is indeed de- 
tected in the Chandra Deep Field South 2-Megasecond 
Catalog (Luo et al. 2008), in both the X-ray channels. 
The luminosities are 3.2 x \(y a crg/s and 4.8x 10 43 erg/s 
respectively in the soft and hard band. On the contrary, 
object 19389 is not detected. Object 24279 has a Sersic 
index of intermediate value n~1.6, suggesting that this 
is probably a bulge dominated disk galaxy. Finally, to 
obtain a reasonable fit to object 24626, we had to use a 
combination of two components: one Sersic model with 
a very high Sersic index to reproduce the inner part of 
the galaxy, and a second one with n~l to reproduce the 
outskirts. 

It is also interesting to analyze the residual maps pro- 
duced by subtracting the best-fit model from the orig- 
inal image to investigate the presence of diffuse low- 
surface brightness halos, as recently suggested by some 
(Mancini et al. 2009; Hopkins et al. 2009b). The residual 
maps show some residual noise in the inner 0.5”: this is 
because we are using an averaged PSF, but the actual 
PSF is not homogeneous across the field. However, this 
effect does not produce any systematic bias in the fitting 
procedure, as it can be seen looking at the 4th column of 
Fig. 1. It is evident that in none of our galaxies such a 
halo is observed. The only exception could be the case of 
galaxy 24626, which requires a two-component model to 
reproduce its light profile, an inner bulge-like part and an 
outer disk. Subtracting both of these two components, 
however, leaves a complex residual map, revealing a sub- 
stantially more irregular light distribution compared to 
the other galaxies. Finally, note that both components 
have effective radius r e ~3 kpc, implying that even if 
observed in shallower images that miss the low surface 
brightness component, the effective radius would still be 
found to be ~ 3 kpc. 

To further explore the possible presence of a faint halo 
around these galaxies we have stacked together all of 
them except galaxy 24626 and then fitted the stacked 
image with GALFIT following the same procedure used 
for individual galaxies. The stacked image, which also 
looks very compact, and the GALFIT residual map are 
shown in Figure 2. The residual image, consistent with 
what we have found analyzing the individual galaxies, 
does not show any large-scale diffuse halo (close com- 
panions have not been masked before the stack and have 
not been fitted). For all the galaxies and for the stacked 
image, however, the amount of residual light in the inner 
2” is lower than 2% of the total light in the galaxy within 
the same aperture. 

In Figure 3 we compare the Sersic profiles in the z, Y, 
.J and H bands, sampling the morphology from the rest- 



Fig. 2. — Stack for the 5 compact galaxies (19389, 22704, 23495, 
23555, 24279) galaxies in the sample (left, panel), and residuals after 
subtracting a Sersic model ( right panel). 
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Fig. 3. — For the 5 galaxies in the sample best fitted by a Sersic 
model, we show the Sersic index ( upper panel) and the effective 
radius (bottom panel) as a function of the wavelength. Object 
23495. that is best-fitted by a point like source, is not reported in 
this figure. The typical error bars are shown, and amount to about 
30% on the Sersic indices and to about 20% on sizes. 


frame UV at A ~ 3300 A to the optical one at ~ 6000 A 
(assuming mean redshift < z >= 1.7). Consistent with 
our visual inspection, this shows weak or no evidence of 
morphological k-correction, since the Sersic indices only 
slightly decrease from the z through the H band, im- 
plying only a slighty higher more compact morphology' 
at bluer wavelengths than at redder ones. This is in 
agreement with reports that the morphology of elliptical 
galaxies is largely independent of the wavelength (Pa- 
povich et al. 2003: Cassata et al. 2005). 

Similarly, the effective radius of all the six galaxies 
in our sample does not depend on the wavelength (see 
bottom panel of Fig. 3), with two of them showing no 
dependence of the size on the wavelength within the er- 
rors, and three showing a slight, but statistically signifi- 
cant, decrease with increasing wavelengths, being ~ 40% 
smaller in H band than in the z one. This results in a 
negative color gradient, with the outer part of the galax- 
ies being bluer than the center. This color gradient has 
also been measured by Guo et al. (in preparation) using 
aperture photometry, and they report that SED fitting 
shows that the properties of the stellar populations in 
the outer parts of these passive galaxies are. on average, 










Fig. 4. — The distribution of effective radii versus stellar mass 
for galaxies in our sample (red filled circles). The continuous 
and dashed lines represent the local relation and its scatter, from 
Shen et al. (2003). Galaxy 23495 is shown as an upper limit, 
since its best-fit is a stellar- like profile. Crosses, circles, squares, 
stars, triangles, pentagons and diamonds indicate the early-types 
of Cimatti et al. (2008), Zirm et al. (2007), Daddi et al. (2005), Van 
Dokkum et al. (2008), Trujillo et al. (2006), Mancini et al. (2009) 
and Longhetti et al. (2007). Mass measurements have been ho- 
mogenized to CB07 models and Salpeter IMF. 

~ 0.5 Gyr younger than those in the center. We conclude 
by noting that the mild dependence of the morphology 
of our galaxies on wavelength is in good agreement with 
other studies of early-type galaxies at high-redshift from 
ACS and NICMOS images (McCarthy et al. 2007; Tru- 
jillo et al. 2007). 

4. THE MASS-SIZE RELATION AT 2 ~2 

In Figure 4 we compare the relationship between mass 
and size (in the H band) for the six galaxies in our sample 
to other measures and to local galaxies (see figure cap- 
tion). The mass measurements have been homogenized 
to CB07 and Salpeter IMF, applying typical offsets: we 
used Mcbqi = Mbco3 — 0.2 and Ms a ip = Mchab + 0.25 
(Salimbeni et al. 2009). Four out of our six galaxies are 
located below the local relation for early-type galaxies, 
in agreement with other repoirts that massive galaxies 
at z ~ 2 are significantly more compact than their local 
counterparts. The remaining two (19389, the object with 
a stellar component at the center, and 24626, the large 
elliptical with multiple components) are on the local re- 
lation. We find it intriguing that these two galaxies also 
are the ones with the lowest redshift in our sample. 

Thus, the sample that we have discussed here confirms 
the existence of very compact and dense passive galaxies 
at ~ 2, on average about 3 to 5 times smaller and about 
50-100 times denser than local counterparts, reported by 
other authors from similar but shallower and/or lower 
quality images (Trujillo et al. 2007; Zirm et al. 2007; 
Longhetti et al. 2007; Van Dokkum et al. 2008) or 
at UV rest-frame wavelengths (Daddi et al. 2005; 
Cimatti et al. 2008). It, however, also suggests a diversity 
of morphological properties (size and stellar mass den- 
sity) among these galaxies, with some of them, possibly 
the ones detected at lower redshifts (i.e. 1.3 < z < 1.5) 
being rather similar to those observed in the local uni- 
verse. 


5. CONCLUSIONS 

We have studied the UV and optical rest-frame mor- 
phology of 6 passively evolving galaxies 1.3 < 2 < 2.4, 
four of which have spectroscopically confirmed redshifts. 

The galaxies have on average a relatively smooth and 
regular morphology, which only mildly depends on wave- 
length. Based on visual inspection and on the analysis 
of the azimuthally averaged light profile, which is rela- 
tively well described by a Sersic profile, we classify the 
six galaxies as early types. This is consistent with their 
selection, which is solely based on their broad-band SED, 
suggesting a tight correlation between spectral and mor- 
phological properties in close similarity to what is ob- 
served in the local universe, i.e. the Hubble sequence. 
Two of our galaxies appear to include a spatially unre- 
solved (at the resolution afforded by the combination of 
1 1ST and WFC3/IH) source at the center that could be 
indicative of an AGN. For one galaxy, object 23495, the 
central point source is responsible for the bulk of the ob- 
served luminosity in the H-band, and this object is in fact 
detected with Chandra, having Lx ~ 3.5 x 10 i3 erg/s. 

Overall the Sersic analysis yields similar results regard- 
less of the bandpass indicating that for these galaxies 
the morphological k-correction between the rest-frame 
UV and optical is small. This particular property opens 
the possibility to study the morphology of these galaxies 
using large, deep and well-controlled samples available 
from the optical images of the GOODS and other sim- 
ilar surveys. We plan to report on such studies in a 
forthcoming paper (Cassata et al. in preparation), where 
we compare the distribution of mid-UV morphology of 
~150 passive galaxies at 2 ~ 2 from GOODS (z band) 
with those of local early-type galaxies from the SDSS at 
the same rest-frame wavelength (u band). 

We do not see any evidence of a faint halo surrounding 
the compact core of these galaxies, neither in the individ- 
ual residuals (after subtracting the best-fit Sersic profile), 
nor in stacked images, for five out of six galaxies in our 
sample. One galaxy, object 24626, requires at least two 
components to provide a good fit, one Sersic model with 
high n (spheroid) embedded in a low surface brightness 
disk. Curiously, both components have similar effective 
radius (but different luminosity), both in the z and H 
bands, implying that in shallower surveys that miss the 
extended component the measured effective radius would 
be the same as the one we derived from the very deep 
WFC3 images. 

Four out of six galaxies in our sample lie below the 
local mass-size relation for early-type galaxies, confirm- 
ing previous reports that such galaxies at z~ 2 can be 
~ 3 times more compact at a given mass and thus ~ 50 
times denser than locally. One of our galaxies, how- 
ever (object 19389), which likely has an unresolved cen- 
tral component, is on the boundary of the local rela- 
tion, while another (object 24626) is well within the lo- 
cal relation. Thus, the deep WFC3 images also show, as 
suggested in other previous studies (Cimatti et al. 2008; 
Mancini et al. 2009) that there are early-type galaxies 
at 2 > 1.5 with size and mass similar to those of local 
counterparts. 
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